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Abstract 

We analyze mathematical models of the global human population growth and compare them to 
actual dynamics of the world population and of the world surplus product. We consider a possibility 
that the so-called world's demographic transition is not a dynamic crossover but a phase transition 
that affects all aspects of our life. 
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I. INTRODUCTION 



A. Empirical models 

The human population of the Earth N& attracted much attention after publication of 
the seminal work of Malthus who realized that it should exhibit the unlimited exponential 
growth: 

dNE at m 
— = rN E . (1) 

The fears were partially dispersed by Verhulst who introduced the logistic equation, 

dN ( N\ 

to account for the population dynamics of closed communities. Here, r is the growth rate 
and K is the carrying capacity. This equation accounts fairly well for the growth of small 
communities but it fails to describe the long-time dynamics of the human population of the 
Earth. 

As it was shown in the seminal work of von Foerster, Mora, and Amioi^, the available to 
them data (up to year 1960) could be fairly well described by the empirical dependence 

**« = it^r (3) 

with the parameters: a « 1, C = 1.8 • 10 11 and t c =2026. The corresponding growth rate is: 

^ « (4) 
dt {t c -t) 2 ' 1 ' 

The most striking feature of Eqs. [31T41 is the divergence of Ne and dNz/dt at finite time, 
t c . This indicates that the above equations are inappropriate in the vicinity of t c . Indeed, 
since 1960, the global human population growth deviated from the hyperbolic dependence 
indicated by the Eqs J51T41 In particular, the growth rate -^^f - achieved its maximum value 
of ~ 2.1% in 1962 and then was steadily decreasing. This prompted the search for the func- 
tions that approximate the hyperbolic dependence given by Eqs l3lT4l before year 1960 and 
replace them by smoother dependences after 1960. Several empirically found replacements 
have been suggested, including hypergeometric 2 , overlay of several exponential or logistic 
curves^, hyperexponential^., delayed logistic curves^, and others. The most insightful em- 
pirical approach was suggested by S.P. Kapitza^ who modified the EqJUas follows: 

dNE = C 

dt (t c -ty + T 2 ' y } 



Here, r is a microscopic time scale for which Kapitza took the lifespan of a generation, ~ 
45 years. This modification captures the maximum in the relative growth rate and assumes 
that the human population eventually comes to saturation. The subsequent studies sought 
to justify this empirical approach. 



B. Mathematical models 



1. Models considering the carrying capacity of the Earth 

In order to understand the future trends of the global human population growth, several 
non-empirical mathematical models have been developed. These models aimed to derive 
EqsjnUH from the "first principles". This approach implies that the equations are con- 
sequences of some plausible scenario while the parameters are empirical. Most of these 
models^"^^^ 1 ^ quantified the verbal approach of Boserup, Simon, Jones, etc. who at- 
tributed the accelerating growth of the human population of the Earth, Ne, to positive 
feedback between the population size and the Earth's carrying capacity, Ke- Then, in ad- 
dition to Eqj2] which accounts for the fast growth of the world population, these models 
introduced an additional equation that accounts for the slow dynamics of the population 
growth resulting from the gradual increase of the carrying capacity: 

iJ t = « < 6 » 

The coefficient 7 quantifies the rate with which the human race expands the carrying capacity 
of the Earth. 

In such a way, these models assume two rates of the population growth. The fast rate, 
as derived from Eqj2j is j^dN/dt ~ r; while the slow rate, as derived by the EqEJ is 
j^dK/dt ~ 7iV. As far as r >> the instantaneous value of the population size (Eqj2]) is 
N m K (the index E has been omitted) . Then Eqj6] reduces to 

dN *r2 

This equation describes an autocatalytic process and its solution is given by EqJS] where 
C = 7" 1 , t c = ti + 1/7-ZVj, and tj, Ni are initial conditions. Upon approaching t c , the pop- 
ulation size N increases and the distinction between the slow and fast dynamics eventually 
disappears. In the limit 7^ >> r, the Eqs J2|6l yield exponential population growth, iV oc e rt . 



At present we do not know whether the human population will come to saturation in 
future or will grow continuously, although we want to believe that its growth will be somehow 
limited. The Kapitza's conjecture consists in replacing the EqJ7]by the empirical equation 

dN r 2 . 2l N 

— — = — sin 8 

dt 7 r w 

that describes accelerating and then decelerating population growth, whereas the population 
size is eventually stabilized at = — . The EqJHl describes the dynamic crossover. It 
operates with the minimal number of parameters -r and 7, and is mathematically appealing. 
However, it can't be easily justified and the reasons for the maximum of the growth rate 
and for the stabilization of the population size remain obscure. 

2. Models based on the Gross Domestic Product- GDP 

The carrying capacity, being the important parameter of the demographic models, can not 
be measured directly. The model that does not consider explicitly the carrying capacity was 
developed by M. Kremer-ii who considered the gross domestic product, GDP = N(S + m), 
as the key parameter that determines the slow dynamics of the population growth. Here, 
N is the population size, m is the subsistence level, and S is the surplus product. Kremer 
related the GDP to the level of technological development T, as follows: GDP oc N^ 1 T^ >2 , 
where 0i,02 are the exponents that should be found empirically. In fact, Kremer put onto 
quantitative language the verbal approach that had been developed earlier by Kuznets, 
Boserup, Jones, etc. The key assumption of the Kremer's model is that the growth of GDP 
is spurred by the technology growth. 

The original model developed by Kremer is static while Korotaev, Malkov, and 
Khaltourina 1 - developed a family of dynamic models basing on Kremer's ideas. In the 
framework of these models, the dynamic variable that quantifies the technological develop- 
ment is the surplus product, S. The simplest model considered by Ref.— consists of two 
equations 

dN 

=— = rNS/m (9) 

(JjL 

f t=1 NS (10) 

with two empirical parameters: r is the rate of the population growth, and 7 has now the 
meaning of the average creative ability of a person. The parameter m characterizes the scale 



of the surplus product S (it can be chosen to be equal to subsistence level) and it has been 
introduced here to be consistent with the notation of the EqJTJ 

In such a way, the Eq. Ois the modification of the Eq{TJ while EqfTUl captures the Kremer's 
idea. The relation to Kremer's work is even more evident if we notice that for <pi ~ 1, the 
definition of the technological level by Kremer :T oc (GDP/N) 1 ^ 2 , is closely related to the 
definition of the surplus product: (S + m) oc GDP/N . 

The relation between iV and S can be found by dividing Eqj9] by EqJTQl This yields 
N (x S. In such a way, the Egs MTOl describe the positive feedback between the surplus 
product and the population growth - from the one hand; and the positive feedback between 
the increasing population and the growth of the surplus product, from another hancU^. The 
solution of these equations exhibits finite-time singularity for both N(t) and S(t). 

The field of applicability of Egs MTOl is evident from their very structure: the right side 
looks as if it were the first term of the power series in the small parameter S/m. In other 
words, Egs MTOl assume that S/m « 1, i.e. they should describe the period before the 
year 1870 when S/m ~ 1. One can go beyond this approximation and modify the Egs MTOl 
to extend their applicability range to beyond the year 1870. Indeed, if we assume that the 
surplus product goes to creation of new working places, the relation of carrying capacity to 
surplus product is especially simple: K = GDP '/m = N(S/m + 1). We replace EqMby the 
logistic equation where K is expressed through S and find 



This equation introduces negative feedback between the population growth and the growing 
population. This has some stabilizing effect and consequently, the solution of Eqs. [TTfTOl 
does not diverge. In the long run, the growth rate of iV comes to saturation, while the 
growth rate of S is unlimited. The relation between N and 5* can be found by dividing 
EqiTTjlby EqlUJ This yields = T *f. The solution of this equation is 



where Nq = r/j. 

In what follows we analyze the data on the relative growth rates: = ^ : N and 
r$ = ^j- S and compare them to the prediction of Eqs MTTlfTSl Figure [TJ shows that 
the relation between iV and 5* is superlinear and is fairly well described by Eqfl2l with 




(12) 
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Nq = 1.65 • 10 9 , the agreement breaks only for N > 3 ■ 10 9 (this corresponds to years 1960- 
1970). Figure [2] shows that the human population growth rate linearly increases with N, 
achieves its maximum value of r^} ax = 0.021 at N w 3 ■ 10 9 and then starts to decrease. The 
EqJTOl correctly predicts dynamics of r^ at N < 2 • 10 9 and deviates from the actual data 
at higher N. Figure [3] shows that the World surplus product growth rate linearly increases 
with N, achieves its maximum value of r™ ax = 0.04 at iV « 3 • 10 9 and then starts to 
decrease. The EqJlO] correctly predicts dynamics of rs at N < 2 • 10 9 an deviates from the 
actual data at higher N. We conclude that the Eqs JT0]fTTI extend the range of applicability 
of the Korotaev, Malkov, and Khaltourina model 15 from the year ~ 1870 to the year ~ 
1960. Since no new parameters/ variables have been introduced, this extension belongs to 
the same family of models^. 

Figure H] shows population growth rate r^ and the surplus product growth rate rs versus 
surplus product, S/m. Both r^ and rg grow with S, achieve the maximum at S/m ~ 5 — 6 
(this corresponds to N ~ 3 ■ 10 9 ) and then decrease. The model correctly predicts initial 
increase of r^ and rs with S (that corresponds to S/m < 1) while pronounced deviations 
from the model occur for S/m > 2. 



3. Models accounting for the demographic transition 

The maximum in the population growth rate (FigfS]) is usually associated with the "de- 
mographic transition"— which results from the decrease of mortality rate and the subsequent 
birth rate decrease. Detailed statistical studies indicate that the reason for decreasing pop- 
ulation growth in developed countries nowadays is the declining birth rate^ whereas there 
is a strong anticorrelation between the level of woman education and fertility. Korotaev, 
Khaltourina and Malkov 15 captured this by introducing an additional dynamic variable: the 
fraction of literate population, I. They modified the Malthus equation, EqJHl to account for 
the negative feedback between the population growth and the literacy level: 

^ = rNS(l-l). (13) 

The dynamics of the surplus product remained the same (Eq JTOj) while the dynamics of I 
has been described by the following equation, 

f f =aSl(l-l), (14) 



where a is a new empirical parameter. 

When the initial educational level of the population is low, this model predicts accel- 
erating growth of N, S, and I. Eventually, when I comes to saturation, N also achieves 
saturation while S does not saturate and grows exponentially. Therefore, this extended 
model captures the nonmonotonous dependence of r N on N (FigfSJ) but fails to account for 
the nonmonotonous dependence of r$ on N (Figj3j). 

4- Critical assessment of the above models 

The common feature of all previously discussed models is that they describe the growth 
of the World human population growth, GDP, surplus product, literacy, etc. by using 
ordinary differential equations containing first-order time derivatives. In the framework of 
these models, the nonmonotonous dependence of the growth rate on time (demographic 
transition) results from the dynamic crossover. This means that at all times there are 
several factors affecting population growth and these factors operate simultaneously. At 
small population number, N < 3 ■ 10 9 , one factor wins and the growth rate increases with 
iV; while for high population number, N > 3 ■ 10 9 , another factor wins and the growth rate 
decreases with N. When iV ~ 3 • 10 9 the gradual transition from one regime to another 
occurs. 

Several features in actual data challenge this picture. First, the transition from increasing 
to decreasing trend in r N vs iV dependence is very sharp (FigfS]). Second, it is not clear 
why transitions in and occur simultaneously, in 1960-1970 and at the same value 
of N ~ 3 • 10 9 (FigHJ). Other parameters also undergo especially fast change in the same 
time period- 1960-1970; these include age structure of the population, the level of literacy, 
urbanization 15 , financial indices 20 etc. All these features are more naturally accounted for 
from another perspective- phase transition. While the commonly accepted approaches 15 
focus on time-dependent, dynamic properties of the population growth, the phase transition 
approach focusses on how demographic and economic variables depend on control parameters 
such as population or surplus product. 
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II. DEMOGRAPHIC TRANSITION AS A PHASE TRANSITION 



The notion of phase transition has been developed in the context of condensed matter 
physics. In the system of many interacting particles/agents, when the control parameter 
(temperature, pressure, density, etc.) varies, the system can progress abruptly from the 
disordered phase where the radius of correlations is finite, to the ordered phase, which is 
characterized by the long-range order. This situation can be usually described using the order 
parameter which is zero in the disordered phase and non-zero in the ordered phase. The 
properties of the system as a function of the control parameter are frequently non-analytic 
at the transition point, in particular, the correlation length diverges upon approaching the 
phase transition and becomes infinite in the ordered phase. Divergence of many physical 
properties at the phase transition is related to the divergence of the correlation length. The 
dynamic properties also undergo dramatic changes and the fluctuations grow on the both 
sides of the phase transition^. The difference between the phase transition and crossover 
scenario is in the following: for the former, the ordered state is characterized by an emerging 
new property- the order parameter which was absent in the disordered state; while for the 
latter scenario all factors were present in both states. 

Phase transitions in social systems (such as financial markets, traffic flow, etc.) were 
noticed long ago^^i. The extrapolated divergence of the human population growth rate at 
years 2026-2040 has been also interpreted as some kind of transitions^. However, the 
growth rate divergence appears only in the "mean- field" models. More realistic models lift 
the divergence and yield earlier date for the population growth switch from one regime to 
another. This probably implies that the phase transition has already taken place in 1960- 
1970, rather than to occur in 2026-2040. Then the demographic transition of 1960-1970 is 
not a purely demographic phenomenon but it is one of the signatures of the global phase 
transition that has been affecting all aspects of human life. 

It is instructive to discuss the properties of this transition of 1960-1970 in the context of 
such a generic phase transition as lattice percolation^. Here, the disordered state contains 
disconnected finite-size clusters while the ordered state is characterized by the appearance 
of the infinite cluster that ensures connectivity of the whole system. This analogy prompts 
us to consider globalization as a hallmark of the phase transition of 1960-1970. One of 
the most prominent aspects of globalization is the economic integration. We consider a 



8 



very crude indicator of the economic integration - the growth of the European Union, in 
particular we focus on r\ = Neu /N Europe -the fraction of the European population in the 
states belonging to European Union (or to its predecessors such as Common Market that 
appeared in 1957). Figure shows dynamics of rj. Amusingly, this very simple measure 
of the European integration mimics the World Globalization index as determined by A. 
Dreher— using weighted economic, political and cultural indicators. 

Figure [5] shows that the onset of European integration can be attributed to the same 
period- 1960-1970 when the global demographic transition occurred. Serrano^ came earlier 
to the same conclusion by considering the historical dynamics of bilateral trade balance. 
Dependence r](t) is very similar to behavior of order parameter at the percolation phase 
transition^. 

If we adopt the hypothesis of the humanity as a system of interacting agents that under- 
goes instability/phase transition peaked at 1960-1970 - this raises many interesting questions 
that prompt extensive scientific research. 

1. What is the nature of the phase transition? What is the difference between two 
phases? Probably, in the years preceding 1960-1970, the most part of the surplus 
product eventually went to the population growth. However, after 1960-1970 the 
surplus product was spent more on the increase of the quality of living, i.e. it was 
channelled to the increase of subsistence level. 

2. What is the proper control parameter that drives this phase transition? Is it the total 
human population, or average population density, or surplus product, or something 
else? 

3. What is the order parameter? Globalization index? Another candidate for the order 
parameter could be some measure of information since the appearance of the global 
information network (Refs.- 1 ^) after 1970 was very prominent. 

4. Which parameters diverge upon approaching the transition? What are the critical 
indices? Besides urbanization and literacy that were studied in Ref.— , it would be 
interesting to consider the historical dynamics of the warfare indicators such as the 
weapon range, the power of the explosives, etc. 
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5. How one can define the correlation length? The possible candidates could be the 
average city population, or the average size of a polity 2 ^. 

6. What is the statistics of fluctuations at this transition? It is well-known that the fluc- 
tuations grow upon approaching the phase transition from both sides. The growth of 
fluctuations in the context of global population growth has been already noticed^. In 
this context, it would be especially interesting to consider the timeline of the financial 
crises. 

7. What is the behavior of the dynamic properties of the global human population at the 
transition? The analog of conductivity of the percolating network would be the signal 
propagation rate or the rate of adoption of technological innovations. How did these 
parameters change through historical time? 

III. SPATIALLY-INHOMOGENEOUS AND DISCRETE MODELS 

Our capability to introduce innovations that enlarge carrying capacity of the Earth (au- 
tocatalycity) translates into human population dynamics equations as a positive feedback. 
It is well known that population dynamics that includes positive feedback and diffusion 
leads to strongly spatially-inhomogeneous population pattern 3 ^ (for example, desert oases, 
vegetation patches in arid zones, etc.) and favors agglomeration. Indeed, the increasing 
economic returns and increasing innovation rate arising from the population agglomeration 
in cities is well documented 2 ^. This means that in the context of the human population 
dynamics the spatial-inhomogeneity by itself has autocatalytic properties. Therefore, the 
models of the human population dynamics should properly address the spatial dimension. 

Note also that the previously discussed models were based on ordinary differential equa- 
tions and didn't take into account neither spatial distribution of the population, nor the 
discrete nature of humans. Very often when the continuum equations describing population 
dynamics assume spatially homogeneous population and predict a very slow growth or even 
population extinction; the individuals self-organize in spatio-temporally localized adaptive 
patches which insure their survival and development. In other words, continuum differen- 
tial equations may fail in predicting the population dynamics of the discrete proliferating 
agents 23 . An interesting example of such approach is the recent study of Yaari et al.— of 



10 



the economics development in Poland after 1990. The Ref.— showed that the economics 
growth was led by few singular growth centers these were associated with the University 
centers. Probably, this shows in a different way the ultimate relation between the education 
level/literacy and the human population dynamics-see Ref.— . 

All this calls for the new generation of the models describing the World human population 
growth. These models should be discrete and spatially-dependent. 

IV. PHYSICAL MEANING OF THE PARAMETERS MEANING OF THE DY- 
NAMIC MODELS 

We consider here a different topic that arises in relation to dynamical models of the human 
population growth. The Egs fTOlfTTl contain two empirical parameters: r, 7 that should be 
somehow related to the human nature. The meaning of the parameter r is more or less 
clear- it is the relaxation rate of the population to sudden changes. It is determined by 
the difference in birth rate and mortality and, to the best of our knowledge, does not 
exceed r record = 0.14. Comparison of the growth rate to the models (FigfSJ) yields r = 
0.013 ~ 0.1r recor rf that is quite reasonable. Note, that r can be measured from the transient 
phenomena, for example, how fast the population size recovers from dramatic disaster such 
as WWII. This yields r ~ 0.02 — 0.03 that is comparable to r = 0.013 determined from the 

The meaning of the parameter 7 is more elusive. Kapitza suggests that 7 = 1/rU 2 where 
U ~67,000 is the coherent population unit. This implies a paradoxical conclusion that any 
coherent population unit consisting of ~ 67,000 individuals will develop into civilization 
consisting of billions of individuals. To get more clear insight into this paradox we compare 
the humans and the beavers. Besides their short lifespan (16-20 years), the beavers remind 
humans in several aspects: they are monogamous, they live in colonies, and most important- 
by building the dams they shape the landscape according to their needs, i.e. they can expand 
the carrying capacity. However, the difference between the "civilization" of beavers and 
human civilization is too obvious. Hence, we seek for deeper relation between the parameter 
7 and the human nature. 

Korotaev et al.— related the 7 to the average creativity of a person. To elaborate further 
on this subject we assume the following scenario. Human population N adjusts to the current 
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carrying capacity very fast. Carrying capacity K slowly increases due to technological 
innovations. So far, the spreading of the technological innovations has been the bottleneck 
that determined the dynamics of the carrying capacity growth. Equation O can be recast as 
follows: 

dK 

= 1 Ndt = 7TdN t (15) 

K 

where N t is the total number of humans lived on Earth till time t, and r is the average 

lifetime. The solution of Eq{T5]is K = Kie" frNt where Ki is the carrying capacity at time U 

and Nt is the total number of humans that lived on Earth between ti and t. Then 

1 dK _ _ x dhiK 
7 ~ 7KdN t ~ T dN t 

For the hyperbolic population growth (EqJSJ), the total number of people lived between ti 

and t depends logarithmically on time: 

N total = -C\n^^- (17) 

The technological innovations are created by people and they are accumulating. This 
means that the carrying capacity at time t is the result of activity of all people that lived 
before^. Therefore, the parameter 7 characterizes the average contribution of an individuum 
to the expansion of the carrying capacity of the Earth. This may be interpreted in two ways. 

1 Each individuum contributes to the growth of the Earth carrying capacity in such a 
way that the average personal contribution is AK = jtK. 

2 The carrying capacity increases by giant steps, AK ~ K, due to stepwise devel- 
opment of science^ and to scientific/technological revolutions^. Then 7T is the 
probability /frequency of these technological revolutions. These revolutions are rare 
events that trigger a series of smaller innovations which become embodied long be- 
fore the next revolution occurs. According to this interpretation, the parameter 7T 
is the probability that an inventor or group of inventors makes a major technologi- 
cal/social/administrative breakthrough. According to this scenario, the human pop- 
ulation growth is a series of logistic curves, each corresponding to a technological 
revolution. The quantity iVo ~ I/7T ( see Eq JT2l) has the meaning of the number of 
people lived that ensure one major technological revolution. 

We believe that the second scenario is more adequate. It has several implications: 
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* Log-periodic oscillations around hyperbolic law given by Eqj3] which were noticed 
by several groups^>22 and attributed to cycles, corresponds to a major technological 
revolution. 

* The current demographic transition and deviation from the hyperbolic law appear 
when the technological revolutions occur so frequently that the full potential of the 
preceding revolution has not been fully realized before the next one occurs. 

* While the motivation for technological innovations so far was the drive towards increas- 
ing carrying capacity, now something changed and the stream of innovations results in 
increased quality of living rather than in increasing number of living persons. (This is 
probably equivalent to increasing subsistence level m). Hence, the population growth 
is not so fast. 

* Is it possible that the very small probability 7T ~ 10~ 9 is somehow related to the 
frequency of genetic mutations which is also exceedingly small (10~ 7 — 10 -8 )? 

* The observation that bigger populations develop fast, while isolated continents, 
archipelagos and islands develop slower may be explained quite naturally. This should 
be related to the probability of appearance of rare events and innovators, and to the 
discreteness of the population. 

The link between the above description and that of Ref.- is provided by the discrete 
character of humans. Indeed, to initiate the positive feedback loop given by Eqj6l for the 
initial group of hominids to expand, it should create at least one working place in the lifetime 
of one generation. This brings us to the minimal group size of iVj = ('jt)~ 1 ^ 2 ~ 67,000. 

Another consequence of the approach based on the number of people lived in the certain 
time interval, is the meaning of "historical time". It has been already noticed 7 that with 
respect to the frequency of historical events, the natural time scale is logarithmic rather than 
linear. Since the total number of humans that lived on Earth also depends logarithmically 
on time (Eq JTTj) . then N t seems to be the "internal clock" of humanity. This conjecture 
provides the basis for quantitative comparison of the historical development of different 
isolated communities. According to this interpretation, the internal clock of a community 
is the total number of people that ever lived in this community. 



13 



Acknowledgments 



I am grateful to Sorin Solomon and Andrey Korotaev for encouragement and for stimu- 
lating discussions. 



1 H. von Foerster, P. Mora, L. Amiot, "Doomsday; Friday, 13 November, A.D. 2026", Science 
132 (1960) 1291-1295. 

2 A. A. Koronovskii, "Population dynamics as a process obeying the nonlinear diffusion equation", 
Dokl. Earth Sciences, 372 (4) (2000), 755-758. 

3 S.D. Varfolomeev, K.D. Gur evich, "The hyperexponential growth of the human population on a 
macrohistorical scale", J. Theor. Biol. 212 (2001) 367-372. 

H. Haberl, H.P. Aubauer, Simulation of human population dynamics by a hyperlogistic time- 
delay equation, J. Theor. Biol. 156, (1992) 499-511. 

5 V.I. Yukalov, E.P. Yukalova, D. Sornette, Punctuated evolution due to delayed carrying capacity, 
Physica D 238 (2009), 1752-1767. 

6 R. Hanson, Long-Term Growth as a Sequence of Exponential Modes (1998), 
|http://hanson.gmu.edu/ longgrow.pdf 

7 S.P. Kapitza, "The phenomenological theory of the world population growth", Sovite Physics- 
Uspekhi, 166 (1996) 63-80. 

8 M. Artzrouni, J. Komlos, "Population growth through history and the escape from the Malthusian 
trap: a homeostatic simulation model", Genus 41 (3-4) (1985) 21-39. 

9 J.E. Cohen, "Population growth and Earth's carrying capacity", Science 269 (1995) 341-346. 

J. Komlos, S. Nefedov, "A compact macromodel of Pre- Industrial Population growth", Historical 
methods 35 (2) (2002) 92-94. 

1 M. Kremer, "Population growth and Technological change: One million B.C. to 1990", The 
Quarterly Journal of Economics, 108 (3) (1993) 681-716. 

2 A.V. Podlazov "Theory of the global demographic process", Mathematical modeling of Social 
and Economic dynamics, Ed. by M.G. Dmitriev and A. P. Petrov, pp. 269-272, Moscow: Russian 
State Social University, (2004). 

3 O. Galor, O. Moav," Natural Selection and the origin of economic growth"; "Evolution and 



14 



growth", European Economic Review 45 (2001) 718-729. 

14 J.C. Chesnais, The demographic transition : Stages, patterns and economic implications, Oxford 
Clarendon Press, 1992. 

15 A. Korotaev, A. Malkov, D. Khalturina, Introduction to social macrodynamics, M. KomKniga, 
2006. 

16 For negative S, Eq(10] leads to a paradoxical conclusion. While the Eqj9] predicts population 
decline for negativs S (that is quite understandable), the Eq{TD] predicts that if S turned out 
to be negative, then it will continue to be negative, in other words S = is the unstable point. 
Following this interpretation, the parameter 7 characterizes the creativity of a person aimed 
towards both constructive and destructive goals! In other words, Eq[10] implies that if some part 
of humanity started some destructive activity, it will pursue it until complete self-destruction. 
Our understanding of the the human behavior would have produced the following equation: 
dS/dt = jN(S + m). Here, if S acquires small negative value (for example, as a result of some 
environmental fluctuation) , then S + m > and the crativity of the humanity brings S back to 
be positive. Only for high enough fluctuation, such that S + m < (as if it were a World war), 
the self-destruction realizes. However, the actual data are in better agreement with EqllOi 

17 D. Stauffer, S. Solomon, Physics and Mathematics Applications in Social Science. Encyclopedia 
of Complexity and Systems Science, Springer 2009: 6804-6810. 

18 B.M. Dolgonosov, V.I.Naidenov, An informational framework for human population dynamics, 
Ecological modelling 198, (2006) 375-386. 

19 C. Chase-Dunn, S. Manning, "City systems and world systems: four millenia of city growth and 
decline", (2002) http:/repositories. odlib.org/irows/ccr021. 

20 A. Johansen, D. Sornette, "Finite-time singularity in the dynamics of the world population, 
economic and financial indices", Physica A 294 (2001) 465-502. 

21 EW. Montroll, WW. Badger, Introduction to Quantitative Aspects of Social Phenomena , 
Gordon and Breach Science Pub., 1975. 

22 A. Maddison, Monitoring the World Economy: A millenial perspective, Paris OECD (2001). 

23 N.M. Shnerb, Y. Louzoun, E. Bettelheim, S. Solomon, "The importance of being discrete: Life 
always wins on the surface", PNAS 97 (19): 10322 (2000). 

K. Sznaid- Weron, "Instabilities in population dynamics", Eur. J. B 16 (2000) 183-187. 
25 R. Taagepera, "Expansion and Contraction Patterns of large Polities", International Studies 

15 



Quarterly, 41 (3) (1997) 475-504. 

26 P. Turchin, Does population ecology have general laws? OIKOS 94: 1726, Copenhagen 2001. 

27 G. Modelski, T.C. Devezas, W.R. Thompson, Globalization as evolutionary process: Modeling 
Global Change , Routledge, 2008. 

28 L.M.A.Bettencourt, J.Lobo, D. Helbing, C. Kuhnert, G.B. West, Growth, innovation, scaling 
and the pace of life in cities, PNAS 104 (2007) 7301-7306. 

29 H. E. Stanley, Scaling, universality, and renormalization: Three pillars of modern critical phe- 
nomena, Rev. Mod. Phys. 71, S358 (1999). 

30 M. A. Serrano, Phase transition in the globalization of trade, J. Statistical Mechanics (2007), 
L01002. 

31 A. Dreher, Does globalization affect growth? Evidence from a new index of globalization, Applied 



Economics, 38, 1091-1110, (2006). [http://globalization.kof.ethz.ch/ 



32 S. Klasen, T. Nestman, Population: population density and technological change, J. Popul. Econ. 
19 (2006) 611-626. 

33 J.B.DeLong, Estimating World GDP, One Million B.C. -Present (http:// 
econl61.berkeley.edu/TCEH/1998-Draft/World-GDP/Estimating-World-GDP.html). 

34 T.S. Kuhn, The Structure of Scientific Revolutions, University of Chicago Press, 1962. 

35 V. F. Turchin, The Phenomenon of Science New York: Columbia University Press, 1977. 

36 D. Stauffer, A. Aharoni, Introduction to Percolation Theory, Taylor and Francis, London 1994. 

37 N.M.Shnerb, P.Sarah, H.Lavee, S.Solomon, Reactive glass and vegetation patterns. Phys. Rev. 
Lett. 90, (2003) 0381011. 

38 G. Yaari, A. Nowak, K. Rakocy, and S. Solomon, Microscopic study reveals the singular origins 
of growth, Eur. Phys. J. B 62, 505513 (2008). 



16 



V. FIGURE CAPTIONS 



Figures 



FIG. 1: Relation between the surplus product S/m and the World human population in the same 
year. The red circles show the data taken from the US Census database and Refs . 11 ' 22 . The 
subsistence level is m = 440 USD. The blue line is the prediction of Eq(12]with No = 1.65 • 10 9 . 

FIG. 2: World population growth rate, r/v = ^jr : N. The red circles show the data taken from the 
US Census database and Refs.— . The blue line is the prediction of the Eq[12]with Nq = 1.65 • 10 9 
and r = 0.013. 

FIG. 3: World surplus product growth rate, r$ = 4jr : S. The red circles show the data taken 
from the Ref. 2 -. The subsistence level is m = 440 USD. The data were averaged over the five-year 
interval. The blue line is the prediction of the Eqllll with Nq = 1-65 ■ 10 9 . 

FIG. 4: Comparison of the World population and the Surplus product growth rates. The data 
for rs were averaged over the five-year interval. Note that maximum and maximum rs are 
achieved simultaneously, at the same value of S/m. 

FIG. 5: Dynamics of the European Union growth, r\ = Neu /^Europe (red circles). Here, Neu is the 
population in the states that belong to European Union or to its predecessors (Common Market), 
while NEurope is the total European population. The blue dashed line shows KOF globalization 
index 31 . Note abrupt growth of r\ around 1960 which is followed by slower steady growth afterwards. 
This is a characteristic behavior of the order parameter at phase transitions. 
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